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Abstract 

The Relativistic Continuum Hartree Bogoliubov theory ( RCHB ), which is the 
extension of the Relativistic Mean Field and the Bogoliubov transformation in the 
' coordinate representation, has been used to study tin isotopes. The pairing corre- 

lation is taken into account by a density-dependent force of zero range. RCHB is 
used to describe the even-even tin isotopes all the way from the proton drip line to 
the neutron drip line. The contribution of the continuum which is important for 
nuclei near the drip- line has been taken into account. The theoretical S2n as well 
as the neutron, proton, and matter rms radii are presented and compared with the 
experimental values where they exist. The change of the potential surface with the 
neutron number has been investigated. The diffuseness of the potentials in tin iso- 
topes is analyzed through the spin-orbital splitting in order to provide new way to 
understand the halo phenomena in exotic nuclei. The systematic of the isospin and 



energy dependence of these results are extracted and analyzed. 
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1. Introduction 

Highly unstable nuclei with extreme proton and neutron ratio are now accessible with the 
help of the radioactive nuclear beam facilities. The extreme proton and neutron ratio of these 
nuclei and physics connected with these unprecedented low density matter have attracted 
worldwide attention in nuclear physics as well as other fields such as astrophysics. The recent 
research in this topic includes the detailed structure, its mode of excitation, and the reaction 
mechanism. Using the measured interaction cross sections it is possible to extract the 
nuclear root mean-square (rms) matter radii, the nuclear mass density and the other related 
observations. In these investigations, the observed sudden rise in the measured interaction 
cross-sections in the neutron rich light nuclei has been attributed to the corresponding large 
increase in the nuclear rms matter radii. This fact is an indication of the sudden change in 
the structure of these nuclei due to the addition of the last few neutrons. The celebrated 
example in these investigations is the discovery of the phenomenon of "neutron halo"0, for 
example in ^^Li. The neutron halo is a state in which neutrons spread like a thin mist around 
the nucleus. It is qualitatively associated with the very small separation energy (Sn) of the 
last one or two neutrons. Due to this small value of Sn in such nuclei, the tail of its last 
neutron wavefunction and so also the corresponding density spreads far out from the center 
of the nucleus. 

New theoretical models and techniques are being developed in order to describe unique 
phenomena in exotic nuclei far from the line of stability. Particularly, on the neutron-rich 
side, the property of exotic nuclei involves the weak binding of the outermost neutrons, the 
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coupling of bound state and the continuum, the diffuse particle distribution and isospin 
dependence of the shell structure. Due to their relevance to the r-process in nucleosynthesis, 
nuclei near the drip-line are also very important in nuclear astrophysics. As for the theoreti- 
cal microscopic description of unstable nuclei is concerned, so far rather different techniques 
have been used, e.g., the Relativistic Continuum Hartree-Bogoliubov ( RCHB ) theoryiH^; 
the exact solution of few-body equations treating inert sub-clusters as point particlesi; the 
Skyrme Hartree-Fock-Bogoliubov (HFB) theorylli; and two-frequency shell modeli. 

The few body method and a shell model diagonalization start by fitting the ground state 
properties of the nuclei or neighboring nuclei and aim at predicting the excited state prop- 
erties. Whereas the relativistic and conventional mean field approaches are based on the 
description of global experimental data for nuclei throughout the nuclear chart. The prin- 
cipal goal here is to obtain a fair description of bulk properties, particularly for the ground 
states. The mean field has certain advantages in describing and predicting quantities, such 
as the radii, mass, tunneling probabilities, etc., in stable nuclei and their future extendability 
to exotic nuclei when the continuum is properly taken into account. 

The RCHB theory, which is the extension of the Relativistic Mean Field ( RMF ) and 
the Bogoliubov transformation in the coordinate representation, provides not only a unified 
description of mean field and pairing correlation but also the proper description for the 
continuum and the coupling between the bound state and the continuum@. 

Recently, a fully self-consistent description of the chain of Lithium isotopes ranging from 
®Li to ^^Li has been given by the RCHbI. The halo in ^^Li has been successfully reproduced 
in this self-consistent picture. Excellent agreement with recent experimental data on the 
radii, binding energy and density distribution is obtained. This remarkable success is mainly 
attributed to the proper treatment of the continuum states and the pairing by solving the 
RCHB equations. This procedure helps to incorporate correctly the scattering of Cooper 
pairs to the 2si/2 level, which is more close to the threshold than the 1^5/2 and 1^3/2 levels, 
in the continuum. The contribution from the continuum is proved to be crucial. Therefore 
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the relativistic mean field provides another picture of halo in nuclei: Due to the extreme 
low proton to neutron ratio, the surface of the neutron potentials become highly diffuse. 
Because the orbital with small orbital angular momentum has less centrifugal barrier, the 
particles with weak binding or even in the continuum in these orbital will have more chance 
to tunnel. In order to study the influence of correlation and many-body effects, it would be 
very interesting to find also nuclei with a larger number of neutron distributed in the halo. 
Based on the RCHB, a new phenomenon "Giant Halo" has been predicted in the Zr nuclei 
close to the neutron drip line, which are composed not only by one or two neutrons, as is 
the case in the halos investigated so far in light p-shell nuclei, but which contain up to 6 
neutron^. 

The development of a proton skin as well as neutron skin in Na isotopes has been sys- 
tematically studied with RCHB in Ref.i, where the pairing and blocking effect have been 
treated self-consistently. A Glauber model calculation has been carried out with the den- 
sity obtained from RCHB. A good agreement has been obtained with the measured cross 
sections for ^^C as a target and a rapid increase of the cross sections has been predicted for 
neutron rich Na isotopes beyond ^^Na. After systematic examination of the neutron, proton 
and matter distributions in the Na nuclei from the proton drip-line to the neutron drip-line, 
a relation between the tail part of the density and the shell structure has been found. The 
tail of the matter distribution is not so sensitive to how many particles are filled in a major 
shell. Instead it is very sensitive to whether this shell has been occupied or not. The physics 
behind the skin and halo has been revealed as a spatial demonstration of shell effect: simply 
the extra neutrons are filled in the next shell and sub-shell. This is in agreement with the 
mechanism observed so far in the halo system but more general. 

From the mean field point of view, the properties of all the nucleons in the nuclei are 
determined by the mean potential provided by their interaction with the other nucleons. 
Therefore the study of the isospin dependence of the potential, which become highly diffuse 
near the particle drip line, is crucial to understand unstable nuclei and the study of the 
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surface diffuseness will provide another mean to understand exotic nuclei. In the present 
work, we extend our former investigation and study the surface diffuseness, the isospin 
dependence of the potential, and spin-orbital splitting in tin isotopes. For exotic nuclei, as 
orbitals with small orbital angular momentum near the threshold will extend far out from 
the nuclear center, and thus the pairing correlation can couple the bound states and positive 
energy states. Particularly on the neutron-rich side, as neutron carries no electric charge, the 
neutron drip line is located very far from the valley of /^-stability. It is crucial to understand 
how the isospin dependence of the potential and the diffuseness develops in unstable nuclei. 
The surface diffuseness can be strongly related with the spin-orbital splitting, which could 
be measured experimentally. 

In the present study, isospin dependence of the potential diffuseness, their relation with 
the spin-orbital splitting , and the ground-state properties of drip-line systems are investi- 
gated by means of the self-consistent RCHB approach. The shell structure and its isospin 
dependence in tin isotope is discussed in the canonical basis. The isospin-dependence and 
energy-dependence of the neutron spin-orbit splitting are discussed. An outline of the RCHB 
formalism is briefly reviewed in Sec. II. In Sec. Ill, the isospin-dependence and energy- 
dependence of the neutron spin-orbit splitting, the neutron, proton, and matter density 
distributions are presented for Sn isotopes. A brief summary is given in the last section. 

2. RCHB Theory 

The RCHB theory, which is the extension of the RMF and the Bogoliubov transformation 
in the coordinate representation, is suggested in Ref.i, its detailed formalism and numerical 
solution can be found in Ref.i and the references therein. The basic ansatz of the RMF 
theory is a Lagrangian density whereby nucleons are described as Dirac particles which 
interact via the exchange of various mesons and the photon. The mesons considered are the 
scalar sigma (a), vector omega {uj) and iso- vector vector rho (p). The rho {p) meson provides 
the necessary isospin asymmetry. The scalar sigma meson moves in a self-interacting field 
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having cubic and quadratic terms with strengths g2 and gs respectively. The Lagrangian 
then consists of the free baryon and meson parts and the interaction part with minimal 
coupling, together with the nucleon mass M, and m^, g„, m^, g^, nip, gp the masses and 
coupling constants of the respective mesons and: 

C = - M)V + \d,,ad^a - U{a) - 

+ |m2a;^a;^ - \Rp,R^'' + \mlp^p^ - iF^.F^^^ (1) 
-ga'4'crip - gujipipi) - gpi)pfi)- eijj/iilj. 

The field tensors for the vector mesons are given as: 

^^^u ^ Qt,^ _ gu^ _ gp^pt, X (2) 

For a realistic description of nuclear properties a nonlinear self-coupling for the scalar mesons 
turns crucia 1§: 

U{o) = -m,a + —a' + —a (3) 
The classical variation principle gives the following equations of motion : 

\a■v^yv{f)^^i{M^Vs{f))\^^ = ^ii'i (4) 
for the nucleon spinors and 



(-AfT + U'{a)) = g^ps 
(-A + mDuj^" = gj^'ir) 



A + m^jp^ = gpJ^'{ 



(5) 



^ -A <(f) = ejjiir) 

for the mesons, where 

Vv{r) = g^ip + gppT + ^e{l-T3)4, 

^ (6) 

^s{r) = gacr{r} 



are the vector and scalar potentials respectively and the source terms for the mesons are 



where the summations are over the valence nucleons only. It should be noted that as 
usual, the present approach neglects the contribution of negative energy states, i.e., no- 
sea approximation, which means that the vacuum is not polarized. The coupled equations 
Eq.(^ and Eq.@ are nonlinear quantum field equations, and their exact solutions are 
very complicated. Thus the mean field approximation is generally used: i.e., the meson 
field operators in Eq.(^) are replaced by their expectation values. Then the nucleons move 
independently in the classical meson fields. The coupled equations are self-consistently 
solved by iteration. 

For spherical nuclei, i.e., the systems with rotational symmetry, the potential of the 
nucleon and the sources of meson fields depend only on the radial coordinate r. The spinor 
is characterized by the angular momentum quantum numbers /, j,m, the isospin t = ±- for 
neutron and proton respectively, and the other quantum number i. The Dirac spinor has 
the form: 



r 

F!Hr) , 



Xtit), (8) 



where Yj^{9, 0) are the spinor spherical harmonics and G^/{r) and Fl\r) are the remaining 
radial wave function for upper and lower components. They are normalized according to 

dri\Gnr)\' + \F^\r)\') = 1. (9) 

The radial equation of spinor Eq. (|D can be reduced as : 

e^G'^ir) = (-^ + -)FHr) + (M + Vsir) + VV(r))Gl^'(r) 

dr r (10) 

e^F!^(r) = {+^ + ^)Gnr) - {M + Vs{r) -Vy{r))F!^{r). 
or r 



where 



-(j + 1/2) for J =1 + 1/2 
+ (j + l/2) for J =1-1/2. 

The meson field equations become simply radical Laplace equations of the form: 



(11) 



are the meson masses for (p = a,u, p and for photon ( = ). The source terms are: 
—QaPs — 920'^ {r) — (73(7^ (r) for the a field 

(12) 



s^{r) 



9loPv 

epc{r) 



for the u field 
for the p field 
for the Coulomb field, 



47rrVs(r) = Ef=i(|G.(r)p - |i^,(r)p) 

47rrV.(r) = Eti(|G.(r)|2 + |F,(r)|2) 

47rrV3(r) = E.lidG.Wp + |F,(r)p) - Eti(|G„(r)p + |F„(r 

4vrrVc(r) = E,li(|G,(r)p + |F,(r)p). 



The Laplace equation can in principle be solved by the Green function: 



(r) = / r drGJr,r )sJr ), 

JO 



where for massive fields 



GJr^r') 



1 1 
2m^ rr' ' 



(•g-Ti^k-'-'l _ g-"i0|r-+r'l 



and for Coulomb field 



(13) 



(14) 



(15) 



(16) 



1/r for r > r' 
1/r' for r < r'. 

The Eqs.(|r^) and (|TT|)could be solved self-consistently in the usual RMF approximation. 
For RMF, however, as the classical meson fields are used, the equations of motion for nucleons 
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derived from Eq. (|l]) do not contain pairing interaction. In order to have pairing interaction, 
one has to quantize the meson fields which leads to a Hamiltonian with two-body interaction. 
Following the standard procedure of Bogoliubov transformation, a Dirac Hartree-Bogoliubov 
equation could be derived and then a unified description of the mean field and pairing 
correlation in nuclei could be achieved. For the details, see Ref.@ and the references therein. 
The RHB equations following: 



/ 



^ h-X A ^ 



A -h + X 



where 

h{f, f^) = [a- p + W(r) + /3(M + Vs{r))] 5(f, r") (18) 

is the Dirac Hamiltonian and the Fock term has been neglected as is usually done in RMF. 
The pairing potential is : 

Akk'ir, n = J d?r, J dh[ V,^^,(rr; r"!). (19) 

kk' 

It is obtained from one-meson exchange interaction ^i^i) in the pp-channel and 

the pairing tensor k = V*!!"'" 

Kkk'iry) =< \akak'\ >= Vy(r)*V'^'(f)^ (20) 
The nuclear density is as following: 

p{fy) = j2^'^'i^*^v\n. (21) 

i 

As in Ref.@, 

^kk' kk' ^^^^ pairing potential in Eq. ( [T9| ) is either the density-dependent 

two-body force of zero range with the interaction strength Vq and the nuclear matter density 
Po- 

V{ru r,) = Vo6{r, - r^)^ [1 - cr.a^] (l - ^] (22) 
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or Gogny- type finite range force witli tlie parameter /Xj, Wi, Bi, Hi and Mi {i = 1,2) as the 
finite range part of tlie Gogny force0: 



V{ri,r2) = ^ 



((ri-ra)/^.)' 



{Wi + B.P'' - HiP^ - M.P^P^) 



(23) 



i=l,2 



A Lagrange multiplier A is introduced to fix the particle number for the neutron and proton 
as = Trp„ and Z = Trpp . 

In order to describe both the continuum and the bound states self-consistently, the 
RHB theory must be sovled in coordinate representation, i.e., the. Relativistic Continuum 
Hartree-Bogoliubov ( RCHB ) theoryS. It is then applicable to both exotic nuclei and normal 
nuclei. In Eq. (0), the spectrum of the system is unbound from above and from below 
the Fermi surface, and the eigenstates occur in pairs of opposite energies. When spherical 
symmetry is imposed on the solution of the RCHB equations, the wave function can be 
conveniently written as 



u 



\ r 



la ■ r 



( . Gv\r) \ 



\ r 



cr ■ r 



YUdA)xt{t). 



(24) 



The above equation Eq. ([l7D depends only on the radial coordinates and can be expressed 
as the following integro-differential equation: 



+ -Gu{r) -{E + X- W(r) + Vs{r))Fu{r) + r / r'dr'A{r, r')Fv{r') = 

dr, , r 



dGu{r 
ir 

'El^-lFjj{r) + {E + \~Vv{r)-Vs{r))Gu{r)+r j r'dr' /\{r,r')Gy{r') = 
dr , r 

K 



dGvir 
dr 

^^M!2 _ lFy{r) -{E-X + Vv{r) + Vs{r))Gv{r) + r I r'dr'A{r, r')Gu{r') = 0, 
dr r 



+ -Gv{r) + (E - A + W(r) - Vs(r))Fy(r) + r / r'dr'A{r, r')Fu{r') = 
dr r 



(25) 



where the nucleon mass is included in the scalar potential Vs{r). Eq.(^), in the case of 
5-force Eq.(^2]), is reduced to normal coupled differential equations and can be solved with 
shooting method by Runge-Kutta algorithms. For the case of Gogny force, the coupled 
integro-differential equations are discretized in the space and solved by the finite element 
methods, see Ref.@ Instead of solving Eqs. (|TU|) and (|TT]) self-consistently for the RMF case. 
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now we have to solve Eqs.(p5D and ([TlD self-consistently for the RCHB case. As the calcu- 
lation for Gogny force is very time-consuming, we use them for one nucleus and fixed the 
interaction strength in 5-force Eg. (|2^) . 

3. Results and discussion 

The procedure to solve the RCHB equations is the same as in Ref.i, we solve the RCHB 
equations in a box of the size i? = 25 fm and a step size of 0.1 fm for the parameter set 
NLSnlli self-consistently. The 5-force is used in the pairing channel and its strength is 
properly fixed by the Gogny force as in Ref.i. The number of continuum taken into account 
is decided by a cutoff energy, i.e., only the levels lying within 120 MeV from the Fermi level 
are taken into account. For po "we use the nuclear matter density 0.152 fm~^. The strength 
Vq of Eq. p^) is determined by adjusting the corresponding pairing energy — ^TrAfi; to that 
of the Gogny force D1S0 in ^^^Sn. 

In Fig. 1, we show the two neutron separation energies S'2„ of even tin isotopes as 
a function of the neutron number from the proton drip-line to the neutron drip-line, 
including the experimental data ( solid points ), and RCHB with 5-force ( open circles ). 
The RCHB calculation reproduces the experimental data quite well. They show a strong 
kink at = 50 and A^ = 82. The drip line nucleus is predicted at ^''^Sn in the present 
calculation. The available empirical data of have been well reproduced. 

The magnitude of the kink at A^ = 82 shows that the shell effects in the RMF theory 
are observed to be strong at A^ = 82. Measurement of the mass of ^^'^Sn, in future, would 
clarify the nature of the shell effects at A^ = 82. Apart from the kinks observed at A^ = 50 
and A^ = 82, S'2„ gradually decrease with A^. The S2n near drip line nuclei ^''^Sn decreases 
relatively faster than the case in light nuclei, e.g., Na isotopesil. This is due to the effects 
of centrifugal barrier. In Sn isotopes, the orbital lii3/2 is just above the threshold for nuclei 
near the drip line. Due to its big centrifugal barrier, particles filled in 1^13/2 do not contribute 
to the binding, although the potential become highly diffuse. In light nuclei like Na isotopes, 
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the orbital 2pi/2 and 2^3/2 are just above the threshold. They gain binding energy due to 
the diffuseness of potential with A^. Therefore the total binding energies of the Na isotopes 
remain unchanged with near the drip line. 

The neutron, proton, matter, and charge radii for all even Sn isotopes are given in Table 
1. The neutron and matter radii increase with the neutron number. Also the proton radii 
increases slowly after ^°°Sn. 

Encouraged by the success in describing the 5'2„ for tin isotopes, we proceed to examine 
the spin-orbital splitting for the whole isotopes. As an example, spin-orbital splitting 

'21 + 1 ^ ' 

versus the binding energy: 

21 + 1 ^ ^ 

in ^^°Sn, ^^°Sn, ^^''Sn, ^"^^Sn, ^^°Sn, ^^°Sn, and ^^°Sn are given in Fig. 2 for the neutron 
spin-orbital partners (1^3/2, 14/2), (lfi'7/2, 15-9/2), (1^11/2,1^13/2), (1^1/2,1^3/2), (1/5/2,1/7/2), 
and {lhg/2, 1/111/2), and the proton spin-orbital partners (1^3/2, 1*^5/2), and (I/5/2, 1/7/2)- It 
is very interesting to see that the spin-orbital splitting for the neutron and proton is very 
close to each other, at least for (1^3/2, lc?5/2), and (I/5/2, 1/7/2) cases. The splitting decreases 
monotonically from the proton drip line to the neutron drip line for all the partners. To see 
the underlying reason of this behavior, it is very helpful to examine the origin of the spin- 
orbital splitting in the Dirac equation. For the Dirac nucleon moving in a scalar and vector 
potentials, its equation of motion could be de-coupled and reduced for the upper component 
and the lower component, respectively. If it is reduced in the lower component, it will be 
related with another interesting topic - the pseudo-spin symmetry discussed in Ref.lll and 
the references therein. For the moment, we are interested in the spin-orbital splitting, the 
Dirac equation can be reduced for the upper component as following: 

1 d{2M -Vv + Vs) d ^^ij 



'dr^ E + 2M-Vv + Vs dr dr 
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_ M1 + k) 1 K d{2M-Vv + Vs) ij 

^ r2 E + 2M-Vv + Vsr dr ^ 

^-{E + 2M -Vv + Vs) {E-Vv- Vs)G^i (r) (28) 
where 

-/ - 1, + 1/2 

« = < ' (29) 

I, J = l- 1/2 



The spin-orbital sphtting is due to the corresponding spin-orbital potential 

1 K d{2M -Vv + Vs) 



(30) 



E + 2M -Vv + Vsr dr 

with some proper normalization factor. It is seen that the spin-orbital splitting in RMF 

is energy dependent and they depends on the derivative of the potential 2M — Vy + Vs as 

well as the particle distribution. Therefore we introduce the so-called spin-orbital potential: 

Vis = ^ ^(^^- ^^ + ^^) . The Vis for i^^Sn, ^^Sn, and i^°Sn, are given in Fig.3. The Vj, 
r dr 

for both proton and neutron are almost the same, as the potential Vy — Vs is a, big quantity 
( ~ 700 MeV ), the isospin dependence in the spin-orbital potential could be neglected. 
Therefore the proton and neutron Vis ^-re the same in the present model. That is the reason 
why the spin-orbital splitting for the neutron and proton is very close to each other in Fig. 
2. From ^^°Sn to ^''°Sn, the amplitude of Vis decreases monotonically due to the surface 
diffuseness. 

For the decline of the spin-orbital splitting, it come from the diffuseness of the potential 
or the outwards tendency of the potential. The diffuseness of the neutron potentials Vy + Vs 
are given in Fig. 4 for ^^oSn, ^^ogn, la^Sn, ^^^Sn, ^^osn, leogn, and i^°Sn. It is seen that the 
depth of the potential decreases monotonically from the proton drip line to the neutron drip 
line and the surface of the potential moves outwards. The inserted figures gives the radii 
Rq at which Vy + Vs — —10 MeV as a function of the mass number. The proton potentials 
Vy + Vs for ii0Sn,i20Sn, i^ogn, i^ogn/sogj^^ leog^^^ no^^ ^^^^^ pig 5 ^j^^ the 

increase of the neutron number, the proton potentials are pushed towards outside as well 
due to the proton-neutron interaction. 
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Similar neutron potentials W - Vs for ^oosn, ^lOSn/^oSn, i^ogn, i^osn/sog^^ leog^^^ 
^^°Sn are given in Fig. 6 and its inserted figures gives the radii Rq at which Vy — Vs = 100 
MeV as a function of the mass number. As seen in the above equations, the spin-orbital 
splitting is related with the derivative of the potential Vy — Vs . The surface diffuseness 
happens for both the vector and scalar potential, it occurs in Vy — Vs and Vy + Vs. 

So far we have started the isospin dependence and showed how the spin-orbital split- 
ting is reduced due to the diffuseness of the potential. Another feature of the spin-orbital 
splitting is the energy-dependence given in Eg . (^8|) . The single particle levels in the canon- 
ical basis for the neutron in ^™Sn are given in Fig. 7. The neutron potentials Vy + Vs 
is represented by the solid line and the Fermi level is represented by a dashed-line. The 
depth of the neutron potential in ^''°Sn is about -60 MeV and the single particle lev- 
els begins from -55 MeV for lsi/2. The spin-orbital partners appear in the order of the 
energy (1^3/2,1^1/2), (14/2,14/2), (1/7/2,1/5/2), (2^3/2,2^1/2), (15-9/2, 15-7/2), (24/2,24/2), 
(1/111/2,1/19/2), (2/7/2,2/5/2), (3p3/2,3pi/2), and (1^13/2, 1^11/2), etc. The larger spin-orbital 
splitting for the same / occurs when the particle is located in the middle of the potential. As 
it can be seen, the spin-orbital splitting in (2^3/2, 2^1/2) partners is twice as large as that in 
(lp3/2, lpi/2) and thrice as large as that in (8^3/2, 3j9i/2). Their energy dependence is clearly 
seen. In Fig. 8, the spin-orbit splitting in Eq.(p6l) versus the binding energy in Eq . (|27|) : in 
™Sn are also given in Fig.7 for (^3/2,^1/2), (4/2,4/2), (/7/2,/5/2), (59/2, 57/2), (/i9/2, /iii/2), 
and (^11/2, ^13/2) partners. It gives a strong energy dependence of the spin-orbital sphtting. 
Similar pattern is also appeared in other tin isotopes, but here we just choose ^^'^Sn as an 
example. 

As mentioned above, the spin-orbital splitting is mainly due to the corresponding spin- 
orbital potential Eq.(^), which decreases with the diffuseness of the potential. However, 
for the single particle levels in the same nuclei, the spin-orbital potential is the same, the 
difference is due to the particle distribution relative to the spin-orbital potential. In Fig. 9, 
the spin-orbital potential Vs multiplied by the density distributions for pi/2 orbital in ^''°Sn 
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is given in the upper panel, and the spin-orbital potential and the probability distribution 
of the p-wave in ^™Sn is given in the lower panel. As seen in the lower panel, the particle 
distributions move outside with the increasing energy. The maximum overlap happens for 
(2p3/2, 2^1/2). This clearly explain the features observed in Fig. 8. Another aspect in Fig. 8 
is that the splitting for n = 1 always increases with /. This is because the centrifugal barrier 
keeps the particle away from center and prevent it from running away so that a big overlap 
always happens for larger /. 

Recent developments in accelerator technology and detection techniques around the 
world have made it possible to produce and study the nuclei far from the stability line 
- so called " EXOTIC NUCLEI" . Experiments of this kind have brought about a new per- 
spective to nuclear structure physics: e.g., the neutron halo and neutron skin as manifested 
in the rapid increase in the measured interaction cross-sections in the neutron-rich light 
nuclei01li. These extremely neutron-rich nuclei and physics connected with the low density 
matter in the tails of the neutron and proton distributions have attracted a lot of attention 
in nuclear physics as well as in other fields such as astrophysics. It is therefore of great 
importance to look into the matter distribution and see how the densities change with the 
proton to neutron ratio in these nuclei. As the density here is obtained from a fully micro- 
scopic and parameter free model, it is well supported by the experimental binding energies. 
We now proceed to examine the density distributions of the whole chain of tin isotopes and 
study the relation between the development of halo and shell effects within the model. The 
density distributions for -"^^^Sn, ^-'^'^Sn,-'^^°Sn, ^'^°Sn, ^^°Sn,-'^^°Sn, ^^°Sn, and -"^"^^Sn are given 
for neutron, proton, and matter density in Figs. 10a, 10b, and 10c, respectively. With the 
increase of A^, the density extends toward outside, only a small increase is observed in the 
center. For the proton densities in Fig. 10b, the density for r < 4 fm decreases dramatically, 
and in reward for this a strong increase in the surface is observed with the increasing of the 
neutron number, because the proton number is fixed. This is quite different from the case 
of Na isotopes, where the proton distribution doesn't vary so much with A^ in the surface 
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are The matter densities for Sn isotopes are given in Fig. 10c. Although neutron density 
increases in the center, it is slower than the decrease in proton. Therefore the total matter 
density decreases in the center. Due to the proton and neutron interaction, both proton and 
neutron densities increase in the surface. 

Fig. 11 is the same as Fig. 10, but in logarithm scale in order to see the detailed 
information in the tail. With the increase of A^, the neutron density extends towards outside. 
The increase of the neutron density beyond r = 6 fm shows a gradual variation, except after 
^^°Sn where the tail of the neutron density does not change. The isotopes ^^°~^^°Sn have 
almost the same tail in the density. For the proton, the tail of the proton density in ^°°Sn 
and ^^°Sn shows difference from the isotopes with > 60, which clearly shows that the tail 
is a asymptotic behavior mainly determined by the last filled few particles. For the tail of 
the matter distribution in neutron-rich nuclei, they are mainly determined by the neutron 
density. 

4. Conclusions 

The Relativistic Continuum Hartree Bogoliubov ( RCHB ) theory, which is the extension 
of the Relativistic Mean Field and the Bogoliubov transformation in the coordinate repre- 
sentation, has been used to study the tin isotopes. The pairing correlation is taken into 
account by a density-dependent force of zero range. RCHB is used to describe the even-even 
tin isotopes all the way from the proton drip line to the neutron drip line. The theoretical 
5'2ri, as well as the neutron, proton, and matter rms radii are presented and compared with 
the experimental values where they exist. The change of the potential surface with the 
neutron number has been investigated. The diffuseness of the potentials in tin isotopes is 
analyzed through the spin-orbital splitting in order to provide new way to understand the 
halo phenomena in exotic nuclei. Summarizing the present investigation We can conclude: 

1. The spin-orbital splitting decreases monotonically from the proton drip line to the 
neutron drip line for all / partners. The spin-orbital splitting for the neutron and 
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proton is very close to each other in RCHB, as the potential Vy — V5 is a big quantity 
( ~ 700 MeV ) and its isospin dependence in the spin-orbital potential could be 
neglected. Therefore the proton and neutron Vis are the same in the present model. 

2. For the decreasing of the spin-orbital splitting with N, it come from the diffuseness 
of the potential or the outwards push of the potential. The depth of the potential 
decreases monotonically from the proton drip line to the neutron drip line and the 
surface of the potential moves outwards. With the neutron number A^, the proton 
potentials are pushed towards outside as well due to the proton-neutron interaction. 

3. A strong energy dependence of the spin-orbital splitting appears in tin isotopes. It 
is due to the change of wave function. The maximum overlap between the particle 
distributions and the spin-orbital potential happens when the single particle energy is 
in the middle of the potential. This energy dependent behavior is stronger for lower 
I orbits than the higher I orbits, as the centrifugal barrier in higher I orbits keeps the 
particle away from center and prevent it from running away so that there is always a 
reasonable overlap with spin-orbital potential. 

4. With increase of A^, the neutron density extends towards outside. The increasing of 
the neutron density beyond r = 6 fm shows a gradual variation with N. The neutron 
tail reaches saturation after ^^'^Sn. For the proton, the tail of the proton density in 
^°°Sn and ^^'^Sn shows difference from the isotopes with > 60, which clearly shows 
that the tail is a asymptotic behavior mainly due to the last filled few particles. For 
the tail of the matter distribution in neutron- rich nuclei, they are mainly determined 
by the neutron density. 

5. With the increase of N, the density extends toward outside, only a small increase is 
observed in the center. For the proton, the density for r < 4 fm decreases dramatically, 
and in return for this a strong increase in the surface is observed with the increasing of 
the neutron number, because the proton number is fixed. Therefore, a slight decrease 
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in the center of the total matter density is observed. 
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TABLES 

Table 1. The neutron, proton, matter, and charge radii for Sn isotopes 



A 










A 


rn 


rp 


I'm 


rc 


96 


4.254 


4.393 


4.327 


4.465 


140 


5.154 


4.712 


5.000 


4.779 


98 


4.288 


4.393 


4.342 


4.465 


142 


5.199 


4.730 


5.039 


4.797 


100 


4.320 


4.393 


4.357 


4.465 


144 


5.243 


4.748 


5.077 


4.815 


102 


4.373 


4AU 


4.393 


4.486 


146 


5.287 


4.765 


5.114 


4.831 


104 


4.422 


4.433 


4.428 


4.505 


148 


5.331 


4.781 


5.152 


4.847 


106 


4.469 


4.452 


4.461 


4.523 


150 


5.375 


4.797 


5.189 


4.863 


108 


4.515 


4.470 


4.494 


4.541 


152 


5.419 


4.812 


5.227 


4.878 


110 


4.556 


4.487 


4.525 


4.558 


154 


5.463 


4.826 


5.264 


4.892 


112 


4.601 


4.502 


4.557 


4.573 


156 


5.506 


4.839 


5.301 


4.905 


114 


4.643 


4.517 


4.588 


4.588 


158 


5.549 


4.853 


5.338 


4.918 


116 


4.684 


4.532 


4.619 


4.602 


160 


5.590 


4.866 


5.374 


4.931 


118 


4.721 


4.545 


4.647 


4.615 


162 


5.628 


4.879 


5.408 


4.944 


120 


4.760 


4.559 


4.677 


4.629 


164 


5.661 


4.894 


5.439 


4.959 


122 


4.797 


4.572 


4.706 


4.642 


166 


5.690 


4.910 


5.467 


4.974 


124 


4.833 


4.586 


4.735 


4.655 


168 


5.715 


4.926 


5.492 


4.991 


126 


4.868 


4.599 


4.763 


4.668 


170 


5.739 


4.944 


5.517 


5.008 


128 


4.902 


4.611 


4.790 


4.680 


172 


5.761 


4.961 


5.541 


5.025 


130 


4.934 


4.624 


4.817 


4.692 


174 


5.783 


4.979 


5.564 


5.043 


132 


4.964 


4.636 


4.842 


4.704 












134 


5.013 


4.656 


4.883 


4.724 












136 


5.061 


4.675 


4.923 


4.743 












138 


5.108 


4.694 


4.962 


4.761 
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Figure Captions 



Fig. 1 Two-neutron separation energies of even Sn isotopes as a function of A^, including 
the experimental data ( solid points ) and the RCHB calculation with 5-force ( open 
circles ). 

Fig. 2 The neutron spin-orbit sphtting Eis — '^~'~^2^^ _|_ ^'-^"'"^"^'^^ versus the mass number 
A in tin isotopes for neutron (1^3/2, 14/2), (l5'7/2, l5'9/2), (l«ii/2, 1^13/2), (lpi/2, IP3/2), 
(1/5/2,1/7/2), and (1/19/2,1/111/2) orbital and proton (14/2, 1^5/2), and (1/5/2,1/7/2) 
orbital, respectively. 

Fig. 3 The derivative of the neutron potentials Vv{r) - Vs{r) for ^^°Sn, ^^°Sn, and ^'^°Sn. 

Fig. 4 The neutron potentials Vv{r) + Vs{r) for ^^oSn, iioSn,i20Sn, ^^osn, i^ogn/sogn^ im^^^ 
and ^™Sn. In order to examine the surface diffuseness more clearly, the radii Rq at 
which Vv{Rq) + Vs{Rq) = — 10 MeV has been given as an inserted figure. 

Fig. 5 The proton potentials Vv{r) + Vs{r) for "°Sn,^20Sj^^ isog^^^ i^ogj^^isogj^^ leog^^^ 

^^°Sn. In order to examine the surface diffuseness more clearly, the radii Rq at which 
Vv{Rq) + Vs{Rq) — —10 MeV has been given as an inserted figure. 

Fig. 6 The neutron potentials Vv{r)-Vs{r) for lo^Sn, ii0Sn,i20gn, isog^^ i^og^ isog^^ leog^^ 
and ^™Sn. In order to examine the surface diffuseness more clearly, the radii R^ at 
which Vv{Rq) — Vs{Rq) = 100 MeV has been given as an inserted figure. 

Fig. 7 The single particle levels in the canonical basis for the neutron in ^^''Sn. The neu- 
tron potentials Vv{r) -\- Vs{r) is represented by the sohd line and the Fermi level is 
represented by a dashed-line. 

Fig. 8 The energy dependence of neutron spin-orbital sphtting in ^^°Sn, 
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Fig. 9 The upper panel is the spin-orbital potential Vis multiphed by the density distribu- 
tions of the pi/2 orbital. The lower panel is the spin-orbital potential and the density 
distributions for pi/2 orbital in ^''°Sn. 

Fig. 10 The neutroni (a) , proton (b), and matter (c) density distributions in Sn isotopes. 

Fig. 11 The same as Fig. 10, but in logarithm scale. 
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